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INTRODUCTION 

^  At  present  ultrasonic  techniques,  for  characterizing  defects  are  being  explored  vigor¬ 
ously  at  various  research  centers  around  the  world.  Currently,  the  ultrasonic  approach  to 
defect  sizing  and  characterization  is  beginning  to  find  great  favor  for  in-  or  out-of-service 
use.  Fer°a  recent  review  of  practical  ultrasonic  nondestructive  evaluation  the  rcade^is 
referred^  Silk-  (ifr'  Because  of  recent  advances  in  experimental  ultrasonic  technology, 
increasing  demands  are  being  put  on  quantitative  theoretical  modeling  of  scattering  of 
ultrasonic  (elastic)  waves  by  cracks,  cavities  and  other  material  inhomogeneities  in  an 
elastic  medium.  With  a  view  to  obtain  detailed  information  about  signals  scattered  from 
complex  defects  various  analytical  and  numerical  techniques  have  emerged.  In  this  article 
we  present  a  combined  finite  element  and  eigenfunction  expansion  technique  for  solving 
scattering  problems  involving  multiple  scattering  as  well  as  complex  geometries.  Numeri¬ 
cal  results  obtained  by  this  method  are  compared  with  recent  experimental  results.  In  the 
following  we  first  summarize  the  theoretical  treatment.  Then  we  present  the  finite  element 
and  eigenfunction  expansion  technique  (FEEET).  Numerical  and  experimental  results  axe 
presented  next.  Although  the  results  presented  here  axe  for  homogeneous  materials,  the 
method  can  readily  be  extended  to  composite  materials.  _ _ 

GOVERNING  EQUATIONS  AND  METHODS  OF  SOLUTION 

For  the  development  of  the  main  concepts  behind  solution  techniques  for  scattering  in 
an  elastic  medium  we  will  assume  that  the  medium  is  infinite  in  extent,  is  homogeneous, 
isotropic,  and  linearly  elastic  with  Lame  elastic  constants  A,  n,  and  density,  p.  An  infinite 
solid  right  cylinder  with  material  constants  Ai,  m,  and  pi  is  assumed  to  be  embedded 
in  the  infinite  medium.  For  convenience  it  will  be  further  assumed  that  the  problem  is 
two-dimensional  (plane  strain)  in  planes  perpendicular  to  the  axis  of  the  cylinder.  Let 
xy-plane  be  a  plane  perpendicular  to  the  axis  and  let  the  boundary  of  the  intersection  of 
this  plane  with  the  cylinder  be  C. 

The  equation  of  motion  in  the  elastic  medium  is  given  by 
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where  u(x,  y)e~iut  is  the  displacement  of  a  point  of  the  medium,  r  is  the  associated  stress 
tensor,  and  u  is  the  circular  frequency,  r  is  related  to  the  displacement  gradients  by  the 
relation 

r  s=  AJ(V  •  u)  +  +  uV)  (2) 

where  I  is  the  identity  tensor. 

Now  the  scattering  problem  is  to  find  the  solution  u  of  (1)  such  that 

tt  =  t|l°)  +tt^,(*,y)  outside  C  (3) 

where  u(°)  is  the  incident  field  that  satisfies  (1)  in  the  absence  of  the  inclusion,  and 
u(*)  is  the  scattered  field  that  also  satisfies  (l)  outside  C  and  the  radiation  condition  as 
r  =  y/x2  +  y2  -*  oo.  Furthermore,  if  the  field  inside  C  is  denoted  by  uL,  then  satisfies 
(1)  and  (2)  with  A,  /x,  and  p  replaced  by  At,  /xi,  and  pt,  respectively.  The  boundary 
conditions  on  C  are, 

u  =  on  C 

n  -  t  =  n  •  Ti  on  C  (4) 

Here  it  has  been  assumed  that  the  cylinder  is  perfectly  bonded  with  the  surrounding 
medium  and  n  is  an  outward  unit  normal  to  C.  The  solution,  u(*l,  has  a  surface  integral 
representation  (See,  Varatharajulu  and  Pao  (2)), 

„«»m + /c{«v  •  B'l'Oi c  m 

Here  Is  related  to  the  Green’s  texisor  G  by  the  equation 

=  AJV  •  G  +  m(VG  +  GV)  (6) 

The  primes  indicate  that  the  quantities  are  evaluated  at  points  r'  on  C.  The  Green’s 
tensor  G  satisfies  the  equation 

V  •  f>|0  +  |r')  =  —I6(r  -  r')  (T) 

Note  that  if  r  approaches  C  then  the  left-hand  side  of  (5)  gives  iu(r).  Equation  (5)  is  the 
basis  of  the  various  boundary  integral  approaches  that  have  been  developed  to  solve  scat¬ 
tering  problems  involving  arbitrary  geometry  of  C\  The  common  difficulty  encountered 
when  using  these  methods  is  in  treating  multiple  scatterers,  scatterers  with  sharp  edges, 
and  scatterers  in  inhomogeneous  or  anisotropic  regions.  For  these  reasons  we  have  devel¬ 
oped  a  hybrid  method  that  combines  the  advantages  of  the  finite  element  method  with 
those  of  the  integral  representation  or  of  the  eigenfunction  representation.  The  hybrid 
combined  finite  element  and  eigenfunction  representation  has  been  used  by  us  in  solving 


'4 


several  scattering  problems  that  involve  non-planar  cracks,  multiple  cavities,  and  cracks 
in  a  bounded  medium 


In  this  approach  Equation  (5)  is  solved  by  expanding  u  in  a  set  of  basis  vector 
functions  that  satisfy  (1).  For  the  two-dimensional  problem  under  consideration  these  are 
(Varatharajulu  and  Pao  (2)),  using  plane  polar  coordinates  r,0, 


and 


f  ;  a  =  1 

\  V(4  #„(&!»•) smn0]  ;  <7  =  2 

V  A  [eae\Hn(k2r)cosn0\  ;  <7  =  1 

V  A  [eMeA Hn(k2r)sinn0\  ;  <7  =  2 


(8) 

(9) 


where  e<>  =  1,  and  en  =  2(n  >.0);  Hn(kr)  is  the  Hankel  function  of  the  first  kind. 
The  basis  functions  (8)  and  (9)  are  used  to  represent  the  scattered  field  u(').  On  the 
other  hand,  the  incident  field  u(°)  is  represented  in  terms  of  Reifr*,  where  Re  denotes  the 
real  part  of  Hn ,  which  is  Jn.  The  quantities  &2,  ki  are  the  longitudinal  and  shear  wave 
numbers, 

ki  =  w/ci  ,  kj  =  u>/c2  (10) 

e  i,  c2  being  the  longitudinal  and  shear  wave  speeds  outside  C. 

Thus  one  writes 

»(<»(r)-2>;«  +  B;j!e«i  (U) 

n,<r 

«(,)M  “£[<«+««)  (12) 

n<cr 

Representations  (11)  and  (12)  are  valid  outside  the  scatterer(s) .  So  they  are  valid  also 
outside  an  imaginary  boundary  B  circumscribing  the  scatterer(s) . 

The  field  inside  this  boundary  is  represented  by  finite  elements.  The  two  solutions 
are  then  matched  by  requiring  the  continuity  of  displacements  and  tractions  across  B.  It 
has  been  shown  that  this  finite  element  and  eigenfunction  expansion  technique  (FEEET) 
is  quite  versatile  in  solving  scattering  by  cracks  or  multiple  cavities  in  bounded  media  as 
well  as  structure-medium  interaction  problems  (Shah,  Wong  and  Datta,  (3);  Abduljabbar, 
Datta  and  Shah,  (4);  Shah,  Wong  and  Datta,  (5);  Datta,  Wong  and  Shah,  (6);  Wong, 
Shah  and  Datta,  (7)).  In  the  following  we  describe  the  method  and  present  numerical 
results.  We  also  present  some  experimental  results. 


The  combined  finite  element  representation  inside  B  and  the  integral  representation 
(5)  outside  B  was  used  in  (8)  to  study  dynamic  amplification  of  surface  displacements. 
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HYBRID  FINITE-ELEMENT  AND  EIGENFUNCTION  TECHNIQUE 

For  the  two-dimensional  problem  considered  in  this  paper,  the  solution  for  the 
displacement  v  can  be  divided  into  two  parts:  the  anti-plane  motion  (SH)  in  which 
=  Uy  =  0,  u,  /  0,  and  the  in-plane  motion  in  which  u,  =  0 ,  ux  r  0,  uy  #  0. 
The  SH  wave  scattering  in  semi- infinite  media  and  thick  plates  was  studied  by  Datta, 
Shah  and  Fortunko  (9),  Abduljabbar,  Datta  and  Shah  (4),  and  Abduljabbar  (10). 


Consider  a  surface-breaking  crack  of  length  D  making  an  angle  a  with  the  x-axis  (See, 
Fig.  1).  Note  that  the  crack  may  be  embedded  in  a  material  that  has  different  material 
properties  from  the  surrounding  homogeneous,  isotropic  elastic  medium.  Let  a  plane  SH 
wave  be  given  by 


(13) 


be  incident  from  infinity.  Here  e  =’  k^D,  2  =  x/D,  and  y  =  y/D.  Figure  1  shows  the 
geometry.  The  incident  field  given  by  (13)  will  be  reflected  from  the  free  surface  so  that 
the  total  field  incident  on  the  crack  is 


tt(°)  s  2wocos(eycos'y)et<**m7  = 


(14) 


where  we  have  dropped  the  time  factor  The  scattered  field  will  be  denoted  by 


m 
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The  solution  outside  an  imaginary  surface  B  enclosing  the  crack  can  be  expanded 
as, 

OP 

w**)  =  ^  AnCosnO'  Hnfar')  (15) 

where  r',  9 '  are  defined  as 

r'  =s  (s2  +  y2)  i  ,  r'eosfl'  =  x 

The  surface  B  can  be  chosen  as  a  semi-circle  of  radius  r*  =  Rb.  Note  that  u/1*)  satisfies 
the  stress-free  boundary  condition  on  y  =  0. 

The  region  interior  to  B  (including  B)  will  be  divided  into  finite  elements.  Thus  the 
displacement  in  the  domain  of  an  element  can  be  written  in  the  form 


M 

"i  =££<({, <?)>»' i  (16) 

isl 

where  M  is  the  number  of  nodes  in  the  element,  <•/»•(*  =  l,..,Af)  are  the  nodal  displace¬ 
ments  and  Li  are  the  shape  functions. 

The  variational  formulation  requires  the  minimization  of  the  functional 


f  =  U-  lp,u*  I J  uiTQdxdy  -  ~  +  PiWi )  (17) 

where  Pi(i  =  1, ...,  M)  are  the  nodal  forces  and  the  integral  is  over  the  area  of  the  element. 
U  is  the  strain  energy  given  by 


r.  1  f  ffSujSo  6u)6q,,  , 

°~2 >')  jlTtTt+TyliW 


(18) 


In  the  above  overbar  indicates  complex  conjugation  and  subscript  e  denotes  element 
property.  Minimization  of  /  leads  to  the  equation 


mm  -  (?) 

where  [K]  is  an  MXM  matrix  whose  elements  are  given  by 


„  fl  f *fr ,  dLi  dLj\  ,31*  dL,  dLidL is 

-  L  Lll{^snt) +l llw ^ + w it)+ 


(19) 


(20) 
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For  the  crack  problem  in  the  absence  of  body  forces  under  consideration  the  interior 
nodal  forces  are  zero.  So  assembling  the  element  equations  we  get  the  global  equation 


Sri  Sib 
Sib  Sbb 


where  subscripts  I  and  B  denote  internal  and  boundary  nodes,  respectively.  Here  {pa} 
denotes  the  column  matrix  of  interaction  forces  at  B  between  the  regions  interior  and 
exterior  to  B.  Then  we  require  that 

{wBy  =  {<4°> + «{,*>} 

{Pb}  =  {Pb}  (22) 

where  p*  is  a  weighted  average  traction  derived  from  w  =  u>(°)  +  Details  of  this 
can  be  found  in  Datta  and  Shah  (11).  Using  equation  (20)  in  conjunction  with  (15)  and 
(21),  one  obtains  equations  for  the  determination  of  the  coefficients  An  and  the  nodal 
displacements  and  WB.  Once  A^’s  are  known,  then  the  scattered  displacement  field 
can  be  calculated.  Figure  2  shows  the  amplitude  of  the  scattered  field  (normalized  with 
respect  to  the  amplitude  of  the  incident  field)  on  y  =  0  for  large  fc2x  when  a  =  90°.  The 
incident  field  direction  is  parallel  to  the  x-axis  (7  =  90°).  Also  shown  on  this  figure  are 
experimental  results  obtained  on  a  number  of  calibration  specimens  using  electromagnetic 
acoustic  transducers  (EMAT).  For  details  of  the  experimental  set-up,  reference  may  be 
made  to  Datta,  Shah  and  Fortunko  (9). 

The  method  described  above  can  be  used  also  to  analyze  scattering  in  a  thick  plate. 
Then  the  representation  (15)  for  the  scattered  field  is  replaced  by 


=  £3  BncoaPnyeik*s 


-  53  Cncos0nye~ik**  ,  *{0  (23) 

m0 

where  0n  =  nir/h ,  kn  =  (fcf  —  /3;|)$,  h  being  the  thickness  of  the  plate.  The  fictitious 
boundary  B  is  now  replaced  by  vertical  lines  at  x  =  —xB  and  x  xr  that  join  the 
horizontal  surfaces  y  *  0  and  y  =  h.  Figure  3  shows  backscattered  amplitude  due  to  a 


FIGURE  2.  Scattered  SH  wave  far-field  amplitude  for  a  normal  surface-breaking  crack. 
Grazing  incidence. 


FIGURE  3.  Scattered  SH  wave  amplitude  in  a  thick  plate,  a  is  the  length  of  the  crack, 
h  the  thickness  of  the  plate  and  k,  —  u//c2. 


normal  surface  crack  (a  =  90°).  Also  shown  in  the  figure  are  the  experimentally-measured 
values.  In  the  experiment  both  the  transmitters  and  receivers  were  sensitive  to  n  =  0  and 
1  modes  only. 


Scattering  by  buried  canted  crack  (a  j*  90°)  in  a  thick  plate  has  also  been  analyzed 
using  the  above  technique.  Experiments  have  also  been  performed.  Details  of  the  analyt¬ 
ical  and  experimental  procedures  are  discussed  in  (12).  Figure  4  shows  the  experimental 
configuration  and  the  geometry  of  the  crack.  In  this  case  several  plate  modes  (n=0,  ..., 
4)  were  generated  and  received  by  the  transducers.  Figures  5  and  6  show  the  reflected 
and  transmitted  transient  signals  at  different  receiver  locations  for  the  same  transmitter 
location.  It  is  seen  that  for  certain  locations  the  agreement  between  theory  and  exper¬ 
iment  is  quite  good.  However,  there  are  other  locations  at  which  there  is  significant 
discrepancy.  There  may  be  several  reasons  for  this:  finite  transducer  size,  crack  lengths 
in  the  experimental  configuration  not  being  exactly  the  same  as  that  assumed  for  analy¬ 
sis,  and,  perhaps,  the  most  important  fact  being  that  the  cracks  in  the  experiments  were 
surface-breaking. 

We  now  present  some  results  for  the  plane  strain  case.  In  this  case  the  scattered  field 
u(*)  can  be  represented  in  the  form  (12)  with 

it  =  vtf  + 

«  «  V  A  (**e.)  +  V  A  (*;«.)  .  (24) 
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FIGURE  4.  Experimental  configuration. 

For  convenience,  0*  and  are  now  defined  as  <pnein3,  tjtncin$,  and  the  sum  in  (12) 
is  over  n  only,  but  from  -oo  to  oo.  The  functions  0£,  and  are  such  that 

satisfies  the  stress-free  boundary  condition  on  y  =  0.  Their  expressions  can  be  found  in 
Datta  and  El-Akily  (13).  The  boundary  B  in  Figure  1  can  be  chosen  to  intersect  the  free 
surface  as  shown,  if  the  scatterer  is  very  close  to  y  —  0,  or  it  can  be  wholly  within  the 
half-space  if  the  scatterer  is  embedded  deep  in  the  half-space.  The  latter  was  the  case 
discussed  by  Shah,  Wong  and  Datta  (5),  and  Wong,  Shah  and  Datta  (6).  More  recently, 
scattering  from  surface-breaking  cracks  has  been  discussed  by  Shah,  Chin  and  Datta  (14). 
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FIGURE  5.  Transmitted  pulse  for  a  30°  canted  buried  crack. 

In  Figure  7  we  show  the  results  of  normalized  scattered  vertical  surface  displacements 
due  to  a  buried  circular  cavity.  Results  of  experiments  and  earlier  asymptotic  calculations 
are  also  shown.  Figure  8  shows  the  corresponding  results  for  a  buried  circular  elastic 
inclusion.  It  is  seen  that  in  both  cases  there  is  very  good  agreement  between  theory 
and  experiment.  Finally,  Figures  9-11  show  the  results  for  surface-breaking  straight  and 
branched  cracks  for  an  incident  Rayleigh  wave  propagating  in  the  x-direction. 
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FIGURE  6.  Reflected  pulse  for  a  30°  canted  buried  crack. 


FIGURE  7.  Scattered  vertical  surface  displacement  amplitude  due  to  a  Rayleigh 
propagating  in  the  x-direction.  H  =  depth  of  cavity  center  =  2.6 A,  k^A  =  0.6. 
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FIGURE  8.  Scattered  vertical  surface  displacement  amplitude  due  to  a  Rayleigh 
wave  incident  on  a  buried  inclusion.  Properties  of  inclusion  and  matrix  are, 
respectively,  c',  =  5.18fcm/s,  c'2  =  2.87 km/s,  ct  =  6.4km/s,  c2  =  3.13Jfcm/s. 


FIGURE  9.  Scattered  vertical  surface  displacement  amplitude  for  a  normal 
surface-breaking  planar  crack. 
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FIGURE  11.  Scattered  vertical  surface  displacement  amplitude  due  to  a  normal 
branched  surface-breaking  crack. 


CONCLUSION 

Results  presented  above  show  the  versatility  and  success  of  the  hybrid  modeling 
technique  to  solve  ultrasonic  scattering  by  single  or  multiple  defects  in  unbounded  and 
bounded  media.  Comparison  of  theoretical  and  experimental  results  shows  excellent  agree¬ 
ment.  These  studies  show  that  the  hybrid  modeling  technique  is  capable  of  describing  the 
scattering  process  in  a  variety  of  situations  that  involve  more  than  one  defect  and  non- 
planar  defects. 

The  analysis  that  has  been  presented  above  for  a  homogeneous  isotropic  plate  can  be 
extended  to  layered  and  anisotropic  plates  modeling  laminated  composites.  That  work  is 
in  progress  and  will  be  reported  later. 
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